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Abstract The Wigner distribution function is a quasi¬ 
probability distribution. When properly integrated, it 
provides the correct charge and current densities, but 
it gives negative probabilities in some points and re¬ 
gions of the phase space. Alternatively, the Husimi dis¬ 
tribution function is positive-defined everywhere, but it 
does not provide the correct charge and current densi¬ 
ties. The origin of all these difficulties is the attempt to 
construct a phase space within a quantum theory that 
does not allow well-defined (i.e. simultaneous) values 
of the position and momentum of an electron. In con¬ 
trast, within the (de Broglie-Bohm) Bohmian theory of 
quantum mechanics, an electron has well-defined posi¬ 
tion and momentum. Therefore, such theory provides a 
natural definition of the phase space probability distri¬ 
bution and by construction, it is positive-defined and 
it exactly reproduces the charge and current densities. 
The Bohmian distribution function has many potential¬ 
ities for quantum problems, in general, and for quantum 
transport, in particular, that remains unexplored. 

Keywords Wigner distribution • Husimi distribution • 
Bohmian distribution • Phase space 


1 Introduction 

With our perspective of more than a century, the dif¬ 
ferent attempts to construct phase space distributions 
for quantum phenomena seem a bit surprising. From 
the very beginning, after the work of de Broglie in 1926 
[Ij, the big fathers of the quantum theory knew a for- 
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mulation of quantum phenomena that leads to a true 
phase space distribution. According to de Broglie’s the¬ 
ory, it is possible to define simultaneously a position 
and a momentum for an electron. However, they pre¬ 
ferred conscientiously to renounce this natural path for 
constructing true phase space probability distributions. 
They looked for alternative ways to build a phase space 
distribution within the orthodox interpretation. All or¬ 
thodox attempts have inherent difficulties. For example, 
the Wigner or Husimi distribution functions are quasi¬ 
probability distributions. The prefix quasi accounts for 
such difficulties. In this paper, we argue that de Broglie 
showed the path for a true phase space quantum distri¬ 
bution. 


1.1 Historical debates 

The beginning of the twentieth century brought surpris¬ 
ing non-classical phenomena: Max Planck’s explanation 
of the black body radiation [2], the work of Albert Ein¬ 
stein on the photoelectric effect [3], the Niels Bohr’s 
model to account for the electron orbits around the 
nuclei 0j, etc. In order to provide an explanation for 
the underlying physics of such new phenomena, physi¬ 
cists were forced to abandon classical mechanics to de¬ 
velop novel, abstract and imaginative theories and for¬ 
malisms. 

In 1924, Louis de Broglie suggested in his doctoral 
thesis that matter, apart from its intrinsic particle-like 
behavior, could exhibit also a wave-like one [5]. Three 
years later he proposed an interpretation of quantum 
phenomena based on non-classical trajectories guided 
by a wave field jBJ. This was the origin of the pilot- 
wave formulation of quantum mechanics that we will 
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refer as Bohmian mechanics to account for the follow¬ 
ing work of David Bohm In the Bohmian formu¬ 

lation, an individual quantum system is formed by a 
point-like particle whose local velocity is defined from 
a guiding wave. The concept of an electron trajectory is, 
therefore, intrinsic to the Bohmian theory mm and a 
phase space with well-defined position and momentum 
becomes a natural construction. 

Contemporaneously, Max Born and Werner Heisen¬ 
berg, in the course of their collaboration in Copenhagen 
with Niels Bohr, provided an original explanation of all 
quantum phenomena without the need of trajectories 
mm- This was the origin of the so-called Copenhagen 
interpretation (also known as the orthodox formulation) 
of quantum phenomena. In the Copenhagen interpreta¬ 
tion, an individual quantum system exhibits its wave 
or its particle nature depending on the experimental 
arrangement. However, it is not possible to retain the 
concept of trajectory because the theory itself forbids 
simultaneous definition of position and local velocity 
for an electron. Therefore, strictly speaking, there is no 
phase-space in the orthodox formulation of quantum 
mechanics. 

In spite of these orthodox difficulties, the scientific 
community has done a constant effort to effectively 
develop additional formalisms or theorie^] to define a 
quantum phase space with position x and momentum 
p. The first attempt was due to Hermann Weyl m 
in 1927 and Eugene Wigner [T3] in 1932, who devel¬ 
oped the so called Wigner-Weyl transform between the 
Hilbert space and a new quantum phase space. The in¬ 
dependent works realized by Hip Groenewold |T5^ and 
Joe Moyal [16] lead to a complete formalism where the 
time evolution of a quantum system itself was defined 
in phase space alone (not in the Hilbert space). For a 
proper discussion see also mmmm- The Wigner 
function is a quasi-probability distribution because it 
gives negative values at some regions of the phase space. 
Of course, such negative values do not invalidate the 
physical contents of the Wigner function because it re¬ 
produces the correct charge and current densities when 
properly integrated. In addition, some relevant works 
do also give physical explanation to these negative val¬ 
ues mmmm- Their main argument is that as far 
as we cannot measure these probabilities, they are a 
correct and acceptable mathematical tool. At the end 
of the day, only positive probabilities corresponding 

1 It is obvious that new mathematical elements open the 
path for different conceptual interpretations. The authors’ 
point of view is that Bohmian mechanics and orthodox quan¬ 
tum mechanics are indeed two different theories, while the 
Wigner distribution function is a new formalism belonging to 
the orthodox school. In any case, such distinction is not clear 
at all and it is not relevant for the conclusions of this work. 


to measurable observables will be recorded. In 1940, 
Kodi Husimi[25] developed a non-negative phase-space 
probability function, the so-called Husimi distribution. 
However, in spite of being a non-negative phase-space 
distribution, it is also a quasi-probability distribution 
because it does not reproduce correctly values of the 
charge and current densities when integrated. 

By construction, Bohmian mechanics, orthodox quan¬ 
tum mechanics and Wigner distribution functions give 
exactly the same empirical results for all (non-relativistic) 
quantum mechanical phenomena. Therefore, all of them 
are perfectly valid tools to develop practical solutions 
for quantum problems. Every formalism has pros and 
cons. For example, when dealing with phase space, it is 
unquestionable that the Wigner distribution has had a 
large practical utility in many quantum problems, even 
in our days This special issue devoted to 

the use of the Wigner function in the computation of 
quantum electronics shows its unquestionable success. 
Therefore, the main goal of this paper is not to criticize 
the quasi-probability nature of the Wigner distribution, 
but to point out that Bohmian mechanics provides an 
original quantum theory with well-defined phase space. 

It is possible to construct a true distribution of proba¬ 
bility in phase space. At the present moment, it is not 
evident if such Bohmian distributions can be as useful 
as the Wigner ones, but it seems obvious that this true 
phase space distribution merits to be investigated. 

1.2 Phase space in Classical mechanics 

In classical mechanics, a particle has well-defined po¬ 
sition x[t) and momentum p(t). Then, if we consider 
a large number N of trajectories (for example, by re¬ 
peating the same experiment with different initial con¬ 
ditions), it is possible to define a classical phase space 
distribution F c (x,p,t) by counting the number of tra¬ 
jectories at each point {x,p} of the phase space: 

1 N 

F c {x,p,t) = — ^ S(x — Xi(t))6(p — Pi(t)) ( 1 ) 

i —1 

where Xiit) and Pi(t) are the actual position and mo¬ 
mentum of the i -th particle at time t with i = 1,..., N. 
The evolution of this classical phase space distribution 
can be found by directly solving the Newton laws of the 
N trajectories or by solving the Boltzmann equation. 

At this point, we want to clarify that the discus¬ 
sion done in this paper will be focused on one-particle 
systems. The function F c (x,p,t ) is constructed by re¬ 
peating the experiment N —> oo times (or by deal¬ 
ing simultaneously with N independent particles). All 
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conclusions elaborated in this paper about one-particle 
phase space probability functions can be straightfor¬ 
wardly generalized for many-particle systems, but with 
a large increment of notation complexity that would 
later become irrelevant for our conclusions. 0 For the 
same reason, we will consider only one spatial degree of 
freedom and one momentum degree of freedom for each 
particle. 

1.3 Phase space in Quantum mechanics 

The proper definition of the phase space is not so ob¬ 
vious in quantum mechanics. First, let us notice that 
there are several different but empirically equivalent 
theories that successfully account for all quantum phe¬ 
nomena. Strictly speaking, there is no phase space in 
the Copenhagen ( orthodox ) quantum theory. The local 
position and local momentum do not exist, simultane¬ 
ously, in that theory. 

Let us start by discussing why there is no phase 
space in the Copenhagen ( orthodox ) quantum theory. 
According to Dirac, the observables are represented by 
operators in Hilbert space [2T]. The commutator of these 
operators specifies which quantities can be known si¬ 
multaneously and which not . In the case of the 
position operator x and the momentum operator p, we 
have: 

[x,p\\ip) = xp\ip) -px\ i/j) = ih\ip) (2) 

Therefore the orthodox theory does not support the 
simultaneous knowledge of local positions and momenta. 
Strictly speaking, there is no phase space in the ortho¬ 
dox interpretation of quantum mechanics. Nevertheless, 
within the orthodox quantum mechanics, we will see 
later that the Wigner formalism defines a phase space 
with a well-defined value for the position and the mo¬ 
mentum of particles, through the use of the Wigner- 
Weyl transform. 

1.4 Definition of a probability distribution 

Let us specify what we understand by a well-defined 
probability distribution in the phase space, say Fq(x,p ), 
for a quantum (or classical) system. One desires that 
this probability distribution fulfils the probabilities ax- 
iom:0 

2 Strickly speaking, the phase space of a system of 2 inter¬ 
acting particles is not {x,p} but {xi, X 2 ,pi,pz}. 

3 In the text, all the integration limits are from —oo to oo 
and thus we will not write them explicitly. 


Fq{x,p)>0, (3) 

J J F Q (x,p)dxdp = 1. (4) 

In addition, its marginal distribution should give the 
usual position or charge probability distributionf0 

Q(x) = J F Q (x,p)dp, (5) 

Expression (|5j) is an important quantity in quantum 

transport because it is related to the charge density, 
which is a very relevant magnitude in any self-consistent 
solution of the electron transport. Another important 
quantity built from this distribution is the current den¬ 
sity, which can be expressed as: 

J(x)= J pF Q (x,p)dp. (6) 

Several attempts have been done trying to provide 

such a phase space quantum distribution Fq (x,p). Here¬ 
after we report the Wigner distribution function [14] 
(which is nothing else but a mathematical Wigner-Weyl 
transform of the density matrix) and the Husimi distri¬ 
bution (a smoothed version of the Wigner distribution 
function). The negative values of the Wigner distribu¬ 
tion function in some regions of the phase space avoids 
the possibility of this Wigner distribution function to 
be a true probability function. The Husimi distribution, 
by construction, has non-negative [25{ values. However, 
it does not reproduce correctly the marginal distribu¬ 
tions ([b]) and ©. 

If one is interested in a true probability distribution 
in phase space, it seems appropriate to use a quan¬ 
tum theory that has a well-defined phase space, i.e. 
a theory that explicitly accounts for well-defined po¬ 
sitions and local momenta fTlfgl . Such theory exist. The 
Bohmian theory briefly mentioned in the first para¬ 
graph is the desired theory. [^] In Bohmian mechanics, 
the electron is defined at any time by a position plus 
a wave-function. The wave-function provides the local 

4 For simplicity, we avoid the explicitly consideration of the 
charge q of an electron in Eq. § and Eq. ©■ We notice also 
that the words charge or current densities can be misleading 
when a wave packet is partially transmitted or reflected, and 
the charge is in fact either fully transmitted or fully reflected, 
not both, when measured. 

5 It must be said that, due to the measurement processes, 
position and momentum are not accessible simultaneously in 
a laboratory in a single experiment as seen in expression 
but the Bohmian theory supports the (ontological) definition 
of both quantities. 
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velocity of the particle that provides the electron trajec¬ 
tory when properly integrated. In summary, Bohmian 
mechanics is a theory which provides a well-defined po¬ 
sition and momentum for a particle. Thus, contrarily 
to the orthodox theory, the Bohmian theory allows the 
existence of a physical and natural phase space. By con¬ 
struction, such Bohmian probability distribution is non¬ 
negative and it satisfies all the probability axioms in 
order to be a correct probability distribution. If the sci¬ 
entific community is interested in describing quantum 
phenomena in phase space, why not using the Bohmian 
distribution function ? 

After this large introduction, in Sec. [2] Wigner, Husimi 
and Bohmian distributions are defined and we compare 
their properties. In Sec. [3] we show numerical results for 
the phase space of the mentioned distributions when an 
electron impinges on a double barrier. We discuss and 
show numerically some disadvantages of the Wigner 
and Husimi distributions, which are not present in the 
Bohmian distribution. Finally, in Sec. [4] conclusions are 
exposed. Several appendixes discuss many technical de¬ 
tails omitted in the main text. 


ensemble in the pure state | ipj). For the sake of sim¬ 
plicity we avoid the explicit time dependence of the 
wave function and Cj. Therefore, the Wigner distribu¬ 
tion function for a mixed state is the following: 

F W (x,p ) ^A x + | )e wv/h dy. (8) 

3 

Hereafter, being the extension to mixed states straight¬ 
forwardly E03 achieved without any modification in the 
conclusions reached in this paper, for simplicity, we will 
consider only pure states. 

The Wigner distribution is a quasi-probability dis¬ 
tribution because it does not satisfy, in general, the 
condition given in Eq. (J3| for a well-defined phase space 
probability distribution. In Appendix [B] we can see why 
this distribution can be negative at some regions of the 
phase space. In addition, in Sec. [3] we report a numeri¬ 
cal example where the probability in the phase space is 
clearly negative. 

We can calculate also Qw{x) and Jw{x) directly 
using Eq. ([5]) and Eq. Q. The results obtained (the 
derivation can be seen at Appendix |A| are: 


2 Quantum phase space probability distribution 

Next, we define and compare the Wigner, Husimi and 
Bohmian phase space probability distributions for quan¬ 
tum systems. 


2.1 Wigner distribution 

One quite common way of describing a quantum me¬ 
chanical system in phase space is by the so-called Wigner 
distribution ( F w ). For a given state \ip), one can con¬ 
struct the density matrix operator p = |^)(i/>| and ex¬ 
press it in the position representation (x\p\x')=(x\il>)(il>\x') 
or in the momentum representation (p\p\p')={p\il>) (if\p'). 
Therefore, somehow, the Wigner distribution can be 
interpreted as an intermediate representation between 
this two and it is given by a Wigner-Weyl transform of 
the density matrix: 


Qw(x) = \if(x)\ 2 


(9) 


Jw{x) = \ip(x)\ 


2 dS(x ) 


dx 


( 10 ) 


being S{x) the angle of the polar representation of 
the wave function ip(x,t) = R(x)exp(iS(x)/fi) (for a 
detailed explanation see Appendix |A|) . The results ([9]) 
and (10) are the ones that one will get by directly us¬ 
ing | if). At this point, because of (|9| and (10), we see 
that F\y{x,p) is a good candidate to study quantum 
transport. However, we will see later that it is “danger¬ 
ous” to take seriously the Wigner (quasi) phase space 
when further developing the basic steps described here. 
For example, when including transitions between the 
phase space points {a :,p} and {x,p'} due to the (Fermi 
Golden rule) scattering. 


F w (x,p) = ^ J if(x+ | W{x- ^)e l ™dy. (7) 

where ip(x ) = ( x\ip }. We notice that the (mathemati¬ 
cal) variable p in the Wigner (quasi) phase space {x,p} 
appears due to a Fourier transform of a type of auto¬ 
correlation function of the wave function. In the case 
of mixed states, the density matrix can be written as 
p = c j IV’jKV’jl where Cj specifies the fraction of the 


2.2 Husimi distribution 

The Husimi distribution ( Fh ) is another possible phase 
space distribution built from the Copenhagen school. In 
this case, it does satisfy the condition © by construc¬ 
tion. For this purpose, we use a set of minimum non- 
orthogonal uncertainty states localized in phase space 
(\q. p)) (3DJ. Using them, the Husimi distribution is the 
following: 
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F H {x,p)= 


1 

2ttH 


(x,p\p\q,p). 


( 11 ) 


If we use that p = we achieve the following 

equation where it is reflected the positiveness of the 
Husimi distribution: 


F H (x,p)=^(x,p\i))(ip\x,p) =^\{x,p\ip)\ 2 . (I 2 ) 

However, it is not also a true probability distribution 
because it does not fulfil the marginal property ([5]). It 
can be also seen that the Husimi distribution is just a 
Gaussian smoothecj^] version of the Wigner distribution 

m- 


If , , -jx-x') 2 -(p-p') 2 2s 2 

Fh(x,p)=— F\v(x ,p )e e h 2 dxdp . (13) 

TThJ 

Next, we will calculate the charge and current den¬ 
sities Qh(x) and Jh{x ) similarly as done for the Wigner 
function (to see the complete derivation see AppendixjA]) : 


QH (*^) 


V2' 


n ts* 


\i’{x') | 5 


(x — x') 2 


dx ', 


(14) 


J H (x) = 


1 


R 2 (x') 


dS(x r ) ^-x'p 
dx' 


dx'. (15) 


\J (27TS 2 ) 

We can clearly see, that these results are the ones 
obtained for the Wigner function, but smoothed by a 


Gaussian function. From Eq. (131 we can understand 


why the Husimi distribution does not accomplish Eq. If. 
The broadening of the probabilities changes the mo¬ 
mentum and position distributions. For these reasons it 
is considered also a quasi-probability distribution. The 
difficulties in properly providing the current and charge 
densities are a dramatic drawback for the correct sim¬ 
ulation of quantum electronic devices with the Husimi 
distribution. 


2.3 Bohmian distribution 

In order to explain the Bohmian phase space distribu¬ 
tion, we briefly explain Bohmian mechanics. It is a the¬ 
ory related to waves and particles. The evolution of the 
wave function is a solution of the typical Schrodinger 
equation: 

= ( 16 ) 

6 We notice that the (mathematical) variables x and p in 
the Husimi (quasi) phase space {x,p} are a smoothed version 
of the x and p variables of the Wigner (quasi) phase space. 


The wave function itself defines a local velocity for 
particles, 


„(*,() = -Im^ 

m ip 


h dS(x i) 
m dx\ 


(17) 


Let us emphasize that particles, here, mean point¬ 
like particles. For a more detailed explanation see m- 
With these two equations, particles can be described 
by trajectories which have, as it has been already com¬ 
mented, a definite position and momentum: 


x(t) = x(0) + f v(x,t)dt. (18) 

Jo 

Once we have well-defined trajectories, we can again 
compute the quantum Bohmian phase space distribu- 
tiorQ similarly to the classical case: 


1 


N 


F B (x,p,t) = lim — 'y'8(x-x i {t))8(p-p i {t)), (19) 
N^-oo iV ' 


i=1 


where N is the number of different trajectories of 
an ensemble of experiments, each experiment has a dif¬ 
ferent initial positiorQ In Eq. (19), Xi(t ) is a position 
of the trajectory in Eq. (18) at time t, while Pi(t) = 
mv(xi(t),t) is the momentum of the particle related to 
the velocity in Eq. ( [TT| with the electron mass to. Let us 
emphasize that, by construction, the phase space distri¬ 
bution constructed with Bohmian mechanics is always 
non-negative. The number of Bohmian trajectories with 
momentum p B at the position x must be positive (or 
zero if there are no particles). 

At this point, we want to emphasize an important 
clarification about the experimental measurement of 
the Bohmian position and momentum. Because of ex¬ 
pression (J2|, in a single experiment, one can measure 
the position or the momentum, not both. However, con- 
trarily to the Copenhagen school, the Bohmian theory 
does not renounce to a well-defined definition of posi¬ 
tion and momentum because of the mentioned experi¬ 
mental limitation. In addition, as we will further discuss 


7 We notice that the variables x and p in the Bohmian 
phase space { x,p } are directly defined in the theory itself. 
They are part of the ontology of Bohmian mechanics. For this 
reason, the Bohmian phase space is a natural space, without 
mathematical tricks. 

8 Let us emphasize that the different Xi(t) and Pi(t) are 

not associated to different particles (as we have said, for sim¬ 
plicity, throughout the paper we only deal with single-particle 
one-degree-of-freedom problems), but to different realizations 
of the same experiment. The probability obtained from the 
wave function has exactly the same (ensemble) mean¬ 

ing. 
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in the conclusions, the Bohmian phase space distribu¬ 
tion is experimentally accessible from a proper treat¬ 
ment of the experimental data obtained from a weak 
measurement of momentum and positions in a large 
ensemble of experiments. 

From the Bohmian distribution, we can calculate 
the charge and current densities. 


Qb(x) = \i>(x)\ 2 , (20) 

J B (x)= \^{x)\ 2dS ^. (21) 

As we can see, these results are exactly the same as 
the ones obtained from the Wigner distribution (and 
different from the ones obtained from the Hussimi dis¬ 
tribution). Therefore, the Bohmian distribution is an 
excellent tool to study quantum electron transport. In 
addition, in the next section we will show how the 
Bohmian phase space distribution does not have the 
problems found for the Wigner and the Husimi distri¬ 
butions. 


3 Numerical example 

According to the conceptual discussions in the previous 
sections, here, we provide numerical examples for the 
three mentioned quantum phase space distributions and 
the related charge density and current density. For sim¬ 
plicity, we consider a simple one-dimensional Gaussian 
wave packet impinging in a symmetric double barrier. 
At the initial time ten the wave function of a Gaussian 
wave packet at the left of the barrier is : 


V’OMo) = ( 


27TOQ 


)\ e ik o(*-xo) ea , p 


{X - X 0 ) 2 
4a 0 


( 22 ) 


where ag = 7.5 nm is the initial spatial variance of 
the wave packet, Xq = 100 nm is the initial central 
position and kg = 0.69nm -1 is the central wave vector. 
In addition, for the Husimi evolution, we used also the 
same dispersion: s = 7.5 nm. 

The time evolution of the initial wave packet is com¬ 
puted by numerically solving the Schrodinger equation 


(16). Then, we compute the three quantum phase space 


distributions at three different times corresponding to 
the initial time tg = 0 ps , the time t\ = 0.09 ps when 
the wave packet is interacting with the barrier and the 
time t 2 = 0.3 ps when the interaction is nearly finished 
and the initial wave packet is clearly split into a trans¬ 
mitted and a reflected components. The information 


corresponding to these three times are plotted in Figs. 
EH and [3] respectively. 

Let us start by comparing the evolutions of the 
Wigner, Husimi and Bohmian distributions in Figs. [!]■ 
[3} It is clearly seen that the Bohmian and Husimi dis¬ 
tributions have non-negative values everywhere at any 
time, satisfying clearly the first probability axiom ([ 3 ]). 
At the initial time, the Wigner distribution is also non¬ 
negative, however, in later times at t\ and £ 2 , negative 
values appears in some regions of the phase space. We 
will further discuss such unphysical feature and their 
consequences later. We also want to emphasize that the 
Bohmian distribution has only one value of the veloc¬ 
ity at each position (this is just a consequence that the 
wave function can take only a single-value at each posi¬ 
tion.). In fact, a realistic example in quantum transport 
must deal with open systems. Then, the pure state has 
to be substituted by a mixed state (a sum of condi¬ 
tional wave functions in the Bohmian language) and 
the Bohmian distribution will provide a distribution of 
velocities (each conditional wave function will have its 
own velocity) at each position in a very natural way. We 
just avoid the consideration of mixed states to simplify 
the present discussion. 

Next, we compare the charge and current densi¬ 
ties calculated using Eq. ([ 5 ]) and Eq. ([6]) for the three 
quantum distributions. Let us emphasize again that, 
as discussed in Sec. [2j the values obtained from the 
Wigner and Bohmian distributions are always exactly 
equal. However, the values of Eq. (14) and Eq. (151 for 
the Husimi distribution does not provide the correct 
charge and current densities obtained from the wave 
function. We clearly see that the modulus squared of 
the wave packet (blue lines) in Figjlbj Fig|2b|and Fig{3b| 
are equivalent to the charge density of the Wigner and 
Bohmian distributions, but not to the Husimi one. 

After confirming, from the numerical simulations, 
the main features that we expect from the distributions 
(i.e. the negative values of the Wigner distribution, the 
mistaken results for the charge and current densities for 
the Husimi distribution and the success in both aspects 
of the Bohmian distribution) we further discuss an im¬ 
portant undesired characteristic of the Wigner distri¬ 
bution. After the interaction with the double barrier, 
say at the time t 2 , the initial wave packets splits 

into a reflected t) part and a transmitted ipri £, t) 
part. 
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(b) Wave packet impinging on a tunneling barrier at to- 




(d) Charge density for the three quantum phase space distributions 
at to - 




(f) Current density for the three quantum phase space distribu¬ 
tions at to- 


Fig. 1: Simulation of the (a) Wigner distribution, (c) Husimi distribution and (e) Bohmian distribution at the initial 
time to- (b) simulation of the wave packet impinging on a double barrier, the simulation parameters are: E = 0.09 eV, 
m* =0.2mo, where too is the free-electron mass, the barrier height is 0.2 eV, the barrier width is 0.8 nm and the well 
depth is 3.2 nm (d) and (f) are the charge and current densities for the three phase space distributions, respectively. 

















Bohmian Function ( a.u.) Husimi Function ( a.u.) Wigner Function ( a.u.) 
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(a) Wigner distribution at t\. (b) Wave packet impinging on a tunneling barrier at t.\. 




(d) Charge density for the three quantum phase space distributions 
at t\. 



(e) Bohmian distribution at t\. 



(f) Current density for the three quantum phase space distribu¬ 
tions at t\. 


Fig. 2: Simulation of the (a) Wigner distribution, (c) Husimi distribution and (e) Bohmian distribution at the time 
t \. (b) simulation of the wave packet impinging on a tunneling barrier with the same parameters as in Fig. [l] (d) and 
(f) are the charge density and current density for the three phase space distributions, respectively. 























Bohmian Function ( a.u.) 
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(b) Wave packet impinging on a tunneling barrier at £ 2 . 




(d) Charge density for the three quantum phase space distributions 
at £ 2 . 




(e) Bohmian distribution at £ 2 . 


(f) Current density for the three quantum phase space distribu¬ 
tions at £ 2 . 


Fig. 3: Simulation of the (a) Wigner distribution, (c) Husimi distribution and (e) Bohmian distribution at the time 
£2 • (b) simulation of the wave packet impinging on a tunneling barrier with the same parameters as in Fig. [l] (d) and 
(f) are the charge density and current density for the three phase space distributions, respectively. 
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As seen in the plot Fig|g of the modulus squared 
of the wave function, the reflected part and transmitted 
part are located at both sides of the barrier (not inside). 

Identically, in the Wigner distribution in Fig |3a| we 
easily recognize the spatial locations of Fw{x,p ) from 
the reflected ipR(x,t) and transmitted ipT(x,t) parts. 
However, in addition, there are large non-zero values 
(negative and positive, because when integrating in this 
region of the phase space the result must be zero) of 
Fw{x,p) in the middle of the barrier, x = 150 nm , at 
places where no probability presence of the electron is 
supposed to be according to Fig|3b| The mathematical 
reason of such spurious result is clear from the Wigner- 
Weyl transform. At time the spatially separated re¬ 
flected if)R(x,t) and transmitted ipT(x,t) parts of the 
wave function i^(x,t) have to be spatially displaced a 
distance y/2 to compute Fw{x,p) from Eq. ([8|. In par¬ 
ticular, for the value of y = 2d, we get i/>r(x + d, t) and 
iPt{x + d,t), plus iPr(x — d, t) and — d,t). Finally, 
according to Eq. ([8]), the Fw{x,p ) at time at position 
xq = 150 nm for any value of the momentum p is: 


X ■ 2 pd 

F w (x c ,p) oc ^r(x c + d)$r( x c - d)e l ~. (23) 

since d is the distance between barrier and the reflected 
part, which we consider equal to the distance between 
the transmitted part of the barrier, the product i/>r(xc+ 
d)ip^{xc~d) in (231 is different from zero at the barrier 
region xc- 

In order to clarify the consequences of this unphys¬ 
ical feature described in Eq. (231 (we know from the 
wave function evolution that there are no electron prob¬ 
ability there!), we remind that by integrating the posi¬ 
tive and negative values of Fw{x,p) around xc we will 
reproduce correctly the charge density at this point. 

There is no charge density at position xc = 150 nm 
in Fig |3d| However, if one tries to gives a physical mean¬ 
ing to Fw(x,p) at these points in Fig 3a (as a true 
physical probability distribution of the electron at the 
phase space), one must be very careful. For example, 
let us imagine that we introduce an ad-hoc scatter¬ 
ing term (due to impurities, for example) in the quan¬ 
tum equation of motion of the Wigner function. If such 
ad-hoc term is introduced as a transition from an old 
phase space point {xctP} towards a new point {xc,p f } 
through the Fermi golden rule probability S p y , we are 
moving electrons from places (for example xc = 150 nm) 
where there are no electrons. The mistake in the transi¬ 
tion of probability from regions without electrons is be¬ 
cause we introduce the scattering mechanism by hand 
as an extra ad-hoc term in the quantum equation of 
motion of the Wigner function. Obviously, this spuri¬ 
ous effect will not be present if the scattering mecha¬ 


nism (with the impurity) is introduced directly in the 
Hamiltonian inside the quantum equation of motion. 
The mentioned undesired features will not occur within 
the Bohmian distribution, because as seen in Fig. [3eJ 
we only see non-zero (positive) probabilities at locations 
where the electron may be reflected or transmitted, but 
not in other regions. 


4 Conclusions 


In this work, we discuss the motivations and interests 
for constructing a true (well-defined) phase space dis¬ 
tribution for quantum systems. We have shown that, 
both, the Wigner and Husimi distributions (related to 
orthodox or Cophengauhen quantum mechanics) do not 
satisfy this purpose completely and, in fact, they are 
called quasi-probability distributions. 


The Wigner distribution, as seen in Eq. (40), car¬ 
ries negative values that cannot be supported by a well- 
defined quantum probability distribution in phase space. 
Of course, when properly integrated, correct probabil¬ 
ity distributions are found out for the charge and cur¬ 
rent densities. However, the ad-hoc manipulation of the 
Wigner distribution function can lead to incorrect re¬ 
sults (we saw for instance in Sec. [3] the non-zero values 
at places where there is no electron and we discuss the 
spurious transitions between regions without electrons 
when introducing scattering by hand as an extra term 
in the quantum equation of motion of the Wigner dis¬ 
tribution). The Husimi distribution (a smoothed ver¬ 
sion of the Wigner distribution), by construction, only 
has positive values. However, we cannot obtain success¬ 
fully the charge and current densities. Thus it is not 
a proper probability distribution for quantum electron 
transport. 


In addition, we have presented the Bohmian dis¬ 
tribution (associated to Bohmian mechanics) that ful¬ 
fills all requirements to be a true (well-defined) prob¬ 
ability distribution. It is a positive probability distri¬ 
bution that exactly reproduces the charge and current 
densities used in the development of quantum electron 
transport simulators. These good properties are just a 
consequence of the Bohmian theory, which allows well- 
defined momentum and position, simultaneously, for 
an electron. The summary of the characteristics of the 
three distributions are seen in Table [lj 
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F w 

Fh 

F b 

Positive distribution 

No 

Yes 

Yes 

Get the exact Q(x) 

Yes 

No 

Yes 

Get the exact J(x) 

Yes 

No 

Yes 


Table 1: Bohmian distribution is the only that satis¬ 
fies all requirements. Wigner and Husimi distributions 
are quasi-probability distributions and do not satisfy 
all requirements. 


Finally, we would like to clarify whether or not the 
Wigner and the Bohmian distributions can be mea¬ 
sured. As we stress in the expression ([2| of the intro¬ 
duction, it is not possible to measure the position and 
the momentum of an electron, simultaneously, in a sin¬ 
gle experiment. From this unquestionable result, the 
Copenhagen school renounces to an (ontological) defi¬ 
nition of the momentum at one particular position (i.e. 
renouncing, in fact, to a proper definition of the phase 
space), while the Bohmian school did not renounce to 
it. Nowadays, with our technological advances in care¬ 
fully measuring (i.e. interacting with) a quantum sys¬ 
tem, the Wigner and Bohmian phase space distribu¬ 
tions are, in fact, associable with experiments in a labo¬ 
ratory. In both cases, the “trick” is repeating the exper¬ 
iment many times (not only one) and treating the data 
of such a large ensemble of experiments in a proper way. 
Regarding the Wigner function, it is “measured” using 
the so-called quantum tomography. For example, one 
can measure the electric field quadratures m in many 
different experiments. Using these data and with the 
help of some mathematical expressions, a phase space 
can be reconstructed through the Wigner function dis¬ 
tribution (with its negative values). See an equivalent 
work in Ref. [ 22] However, one must take into account 
that this (ensemble) phase space is not the true phase 
space of a wave function in a single experiment (which 
does not exist in the orthodox quantum mechanics). 

Identically, Bohmian phase space distribution can 
be measured by an ensemble of repeated experiments, 
through the use of weak measurements [531155] . This 
kind of measurement does not perturb so much the wave 
function as a strong (or projector) measurement, but 
the output result has a large uncertainty [551155] . With a 
proper treatment (the so-called weak value) of the data 
of an ensemble of weak and subsequent strong measure¬ 
ments of the positions and velocities, the Bohmian dis¬ 
tribution can be constructed. In this regard, let us no¬ 
tice that the velocities of photons at particular positions 
have been already measured in a double slit experiment 
[35] . Again, we remark that this result comes from an 
ensemble of experiments. Thus, again, it is just an en¬ 
semble phase space. We cannot extract this information 


from one unique measurement. In any case, contrarily 
to the Copenhagen school, Bohmian mechanics allows 
a well-defined probability distribution in phase space 
(according to the rules ©-©)• To be honest, every 
formalism has pros and cons and Bohmian phase space 
distributions is still in its infancy. In any case, in the 
authors’ opinion, such true phase space distribution has 
an enormous potential to be developed and exploited, 
in general, for study any (non-relativistic) quantum sys¬ 
tem and, in particular, for quantum electron transport. 
Some results justifying the viability of this path can be 
found in the BITLLES simulator developed by Oriols 
et al. [5511571155] . 
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Appendix A Q(x) and J(x) derivations 

In this appendix, we will develop the calculus to obtain 
the charge and current densities (Eqs. ®> and ((fif) for 
each of the three analysed distributions. 


A.l Wigner distribution 

The charge distribution is straiglrtforwadly found out: 

Qw{x,p ) = ^ J dp J i/j(x + |)^*(x - | )e l ^dy = 

= J S(y)ip{x + | )-ip*(x - ^)e l ^dy = \^{x)\ 2 . (24) 

The computation of Jw(x) needs a more detailed 
discussion: 


J w {x) = J pF w (x,p)dp = 

= ^ J p J ip(x + | )ip*(x - ^)e l ^dydp. (25) 

Using the chain rule for a function F(x, y) derivable 
and zero-valued at y —> ±oo, we can use the following 
relation 


J dye lpv/n -^F(x, y) = ^p J dye lpv/h F(x,y) (26) 

and using the polar form of the wave function (ip{x) = 
R{x)e lS ( x l), then rewrite Eq. (25) as: 
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—l 


Jw ^ = 2i,J dp d y e ~ ipv/h - 


R( x + y -)R(x- y -). 


r [s(x + #)-s(x-f)] i ( dS{x + %) _ dS(x-%) \ 
h\ dy dy ) 


dy 


r [s(x+jf)-s(*- D] 

dy V 


x +\) R {x~l 


(27) 


To proceed, let us focus on the following term: 


dS(x + |) dS(x—%) v | _ S'(a; + 1 ^) — S(x + |) 


dy 


dy 


2 '-= lim 

t—^ 0 . 


S(x-y±t)-S{x-l)^ 1. r S(x+§ + §)-S(:r+§) 


=- lim 
2 t-y0 L 


S{x-l)-S{x-l- I) 


i/2 


i/2 

dS{x+l) dS(x- 1)- 




9a; 


9x 


,(28) 


which can be used to rewrite Eq. (271 as 


( 


2' v 2 

1^9S(cc+§) 9S(a;-f) 

9a; dx 

^R { x+ y -)R(x- y -))\=R 2 (x) 


^ ihe m*+ D-^-i)]. 

9S'(a;) 


9a; 


(29) 


Eq. (291 is the expression for the current density 
distribution for the Wigner function. 

In Eq. (29) we have used the following property: 

J dpe~ ipv/h = 2nH6(y), (30) 

and the fact that the second term within the integral 
in Eq. (29) is zero: 

IK+I^-I)) 

_ 1 
~ 2 
= 0. 


y=0 


dR(x + §) y y dR{x — 2 ) 


dx 


•-R(x-±)-R{x + ±y 


dx 


y —o 

(31) 


A.2 Husimi distribution 

The derivation of the Husimi charge and current densi¬ 
ties are quite similar to the derivation realized for the 
Wigner distribution, with the only difference that now, 
the Wigner distributions are smoothed with a Gaussian 
function. In first place, we derive the charge distribu¬ 
tion: 


Qh ^ = Ih JJJl j^'-\)r{x' + \)e iP ndy 

(g— e') 2 (p-p') 2 2. 2 1 1 

e e R 2 d x dp dp = —- — • 



(a —a/) 2 


(32) 


After that, we derive the current density. For this, 


Eqs. (26), (28), (30) and (31) will be used again: 


Jh(x)=^JJJpFh{x', p')e 2 . 2 ’ e 2<t p dx'dp'dp= 

= ^ fffn ■ 


_(p-p') 2 2 s 2 1 1 

e s 2 aa; ap dp =— 


h v / (2tts 2 ) 
V>*(a; / + ^)e l£? ^ dye ^ dx'dp = 


pV(^-|) 


1 


V / (2tts 2 ) 


„, . dS(x') (x-x) 2 . 

,2 ^ v ' e ^ rfa;' 


j R\x') 


dx' 


(33) 


Therefore, we can clearly see that these results for 
the charge and current densities are different from the 
Wigner results. 


A.3 Bohmian distribution 

The charge distribution is the following: 


f 1 x , 

Qb(x, t) = / lim —S26(x-Xi(t))6(p-p t .(t))dp = 

/ Af->oo JV z ' 

J »=1 

1 * 

= lim = \ip(x,t)\ 2 . (34) 

Af->oo IV z ' 
i=l 


The last equality in Eq. (34) is due to the quantum 
equilibrium hypothesis, which states that the charge 
distribution at time t is equal to the modulus squared 
of the wave function (for a more detailed discussion see 

ra¬ 


in order to calculate the current density, we need 
to express the wave function in the polar form (in the 
same way as done for the Wigner distribution) and we 
proceed in the following way: 
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1 X , 

J B {x,t)= lim — Y ^p6(x - Xi(t))6(p - pi(t))dp = 

N—too J\J z ' 

1=1 

= J™, !><*(* - *<(*)) = R 2 ( x )^r- < 35 ) 

Apart from using again the quantum equilibrium 
hypothesis, we have also used that the momentum in 
Bohmian mechanics in a single-particle case is 1-‘11)1 : 


f <j)(2x — r)dr = — / (j)(2x — r')dr' = 

J J OO 

= J (j){r')*cj)[2x — r')dr'. (39) 

Therefore, the sum of both terms are exactly equal 


to F\y{x,p) in Eq. (371. 


1 4>(r) + <j>{ 2x — r)\ 2 dr 


\<p(r)\ 2 dr + 


p(x, t ) = Im 


_ dS(x) 
ip dx 


(36) 


Let us notice that Eq. (36) is just Eq. © multiplied 
by the electron mass m. 


+ J\cp(2x — r)\ 2 dr + 2 


(p(r)cp*(2x — r)dr > 0. (40) 


that we rewrite as: 


Fw(x,p) 


| </>(r) + 4>{2x — r )| 2 dr 


Appendix B Negative values of the Wigner 
distribution 

If we manipulate the Wigner distribution (applying a 
change of variable), we can see that it can be thought 
as a correlation function: 


F w (x,p) = ^ J dp J ip{x + ^)ip*(x - ^ )e l ™dy = 
= 2 f ip(r)ip*{2x — r)e‘Te I_Li; “ dr = 2j ip(r)e lSir ■ 

■ ( ip(2x — r)e _l n Ydr = 2 J (j>(r)<j>*(2x — r)dr. (37) 

The change of variable is the following: x + | = 
r. In addition, we have defined 4>{r) = 'ip(r)e l xr. For 
simplicity we do not indicate the dependence on p. 

With these considerations, we can prove that the 
Wigner function is real and can be negative. We com¬ 
pute the integral of a modulus of a function | <j>(r) + 
4>(2x — r)| 2 which is obviously positive: 


J (4>(r) + <j>(2x — r)) • (4>{r)* + (p*{2x — r))dr = 


= J l < / > ( r ) + 4>{2x - r)\ 2 dr=J (|</>(r)| 2 + \(p(2x - r)\ 2 + 
+ (p{r)(f>*{2x — r) + 4>(r)*4>(2x — r)\dr > 0. (38) 


Next, we observe that the last two terms in the right 
hand side are identical and, by using Eq. (37), equal to 
the Wigner distribution function Fw(x,p). in order to 
see that both terms are identical, we make a change of 
variable r' = 2x — r: 


— J \(p(r)\ 2 dr — J \<j)(2x — r)\ 2 dr 


(41) 


Therefore, the value of the Wigner function in Eq. (41) 
must be real and it can clearly take negative values 
when the first integral is smaller than the sum of the 
other two (for example, when </>(r) = —<j)(2x — r)). In 
this way, it is proved that the Wigner function is a 
quasi-probability distribution, but not a true one. 
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